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CONCENTRIC COPLANAR CAPACITIVE 
SENSOR SYSTEM WITH QUANTITATIVE 
MODEL 

CROSS-REFERENCE TO RELATED 5 

APPLICATIONS 

This application claims priority under 35 U.S.C. §119 to 
provisional application Ser. No. 61/365,601 filed Jul. 19, 
2010, herein incorporated by reference in its entirety. 10 

GRANT REFERENCE 

This invention was made with government support under 
Grant Nos. NNX07AU54A granted by NASA and FA8650- 15 
04-C-5228 granted by U.S. Air Force. The Government has 
certain rights in the invention. 

FIELD OF THE INVENTION 

20 

The present invention relates to nondestructive evaluation. 
More specifically, but not exclusively, the present invention 
relates to nondestructive evaluation of dielectric materials, 
including multi-layered dielectric materials, using a capaci- 
tive sensor. 25 

BACKGROUND OF THE INVENTION 

Dielectric materials play an extensive role in both indus- 
trial applications and scientific research areas. In the modern 30 
integrated circuit industry, as electrical components are min- 
iaturized, there are palpable needs for dielectric measure- 
ments of low-loss thin materials. The use of fine-line signal 
conductors requires thinner, possibly laminated, low-dielec- 
tric constant printed-wiring board materials. On the other 35 
hand, compact antenna arrays require high-dielectric con- 
stant substrates to obtain phase shifts. Moreover, lightweight 
structural composites in air- and space-craft, Kevlar body- 
armor and ceramic -matrix-composites for thermal stability in 
hot engine environments are examples of some of the recently 40 
developed applications of low-conductivity materials. As a 
result of these increased applications of dielectrics, the quan- 
titative dielectric property characterization of these dielectric 
materials becomes markedly important for the process con- 
trol in manufacturing, optimization of electrical apparatus 45 
design and performance, and system monitoring and diagnos- 
tics. 

A number of high frequency nondestructive evaluation 
(NDE) techniques have been developed for dielectric mea- 
surements with their own specific applications [ 1 ] . Transmis- 50 
sion-line techniques are capable of measuring material per- 
mittivity by an open-circuit termination. The material 
properties of the test-piece can be interpreted from the reflec- 
tion coefficient of the system. Open resonators have also been 
used in measuring low -loss materials in the millimeter wave- 55 
length range [2] and a certain open resonator system for 
measuring anisotropic thin films has been developed and is 
able to obtain the material tensor permittivity values [3], 
Measurements using surface electromagnetic waves are quite 
applicable for low-loss dielectric thin films and layered sub- 60 
strates, since they possess a high quality factor and are there- 
fore sensitive to loss [4], Evanescent-field dielectrometry has 
been utilized in diagnosing and monitoring fresco degrada- 
tions resulting from moisture and soluble salts [5]. Besides, 
broadband dielectric measurements (0.01 to 3 GHz) on the 65 
effects of exposure of thick film adhesive-bonded structures 
to moisture have been reported [6], where the data obtained 


2 

are complemented by mechanical testing and failure analysis 
of the bond structure measured as a function of the exposure 
time. However, the focus here is on describing electrostatic 
and low frequency NDE techniques for dielectric measure- 
ments. 

One important and practical field of material dielectric 
property characterization is dielectrometry, which derives the 
complex permittivity of a test -piece from the measured sensor 
capacitance. Interdigital dielectrometry sensors, with 
increased effective length and output capacitance between the 
electrodes because of their interdigital structure, have been 
used for dielectrometry measurements for a long time. An 
excellent review paper on interdigital sensors and transducers 
is [7], in which the physical principles, sensor design and 
fabrication, and relevant applications of interdigital sensors 
are discussed in detail. These interdigital dielectrometry sen- 
sors have been applied in many fields such as material prop- 
erty monitoring, humidity and moisture sensing, electrical 
insulation properties sensing, monitoring of curing pro- 
cesses, chemical sensing, biosensing, and so on. For example, 
using a secant method root-searching routine for parameter 
estimation, interdigital electrode dielectrometry has been 
made capable of measuring the continuum parameters of 
heterogeneous media [8], which include material thickness, 
material permittivity with thickness known, and material sur- 
face conductivity with thickness known. The optimization of 
multi-wavelength interdigital dielectrometry instrumentation 
and algorithms has also been described in [9] . Through varia- 
tion of geometrical design, materials, manufacturing pro- 
cesses, electronic circuitry, and considerations of accumu- 
lated effects of non-ideal geometry of experimental setups, 
improvement of sensor performance can be achieved. Addi- 
tionally, design principles for multichannel fringing electric 
field sensors, especially detailed analysis on how the sensor 
geometry affects the sensor performance and tradeoffs among 
different design objectives, have been carried out [10] pro- 
viding insight into design of capacitive sensors in general. 

Apart from using interdigital dielectrometry sensors, other 
sensor configurations have been used to characterize defects, 
moisture content, temperature, aging status, delamination, 
and other inhomogeneities in dielectric materials. For 
example, rectangular capacitive array sensors have been used 
for the detection of surface and subsurface features in dielec- 
trics and surface features in conductive materials [11]. Cylin- 
drical geometry quasistatic dielectrometry sensors with sig- 
nal interpretation based on semi-analytical models have also 
been developed in recent years to measure the permittivity of 
a dielectric plate [12]. For water intrusion detection in com- 
posite structures, rectangular coplanar capacitance sensors 
with high sensitivity have been developed [13] on the basis 
that the presence of defects, such as water, leads to changes of 
dielectric characteristics in the structure, resulting in varia- 
tions in the sensor measured capacitance. Using a similar 
principle, rectangular coplanar capacitance sensors have been 
applied for damage detection in laminated composite plates 
[14]. Also, the influence of electrode configurations on a 
differential capacitive rain sensor, which consists of a sensi- 
tive capacitor whose capacitance changes in the presence of 
water and an insensitive reference capacitor, have been inves- 
tigated in [15]. Moreover, these capacitance techniques have 
even been employed for the continuous monitoring of the 
thickness of biofilms and tissue cultures [16]. 

Electrical capacitance tomography (ECT) is another 
capacitance measurement technique that is used to image 
cross-sections of industrial processes containing dielectric 
materials [17]. The principle is that through image recon- 
struction for ECT, the test-piece permittivity distribution and 
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therefore the material distribution over its cross-section can 
be determined. Over the past decades, research progress on 
both the hardware design [18, 19] and sensor configuration 
optimization [20] of ECT systems has been made success- 
fully. 

Despite these advances in various capacitance measure- 
ment techniques, problems remain. What is needed is a sensor 
and associated methods and systems which can be used in 
applications, such as, but not limited to quantitative charac- 
terization of material properties of multi-layered structures, 
detection of water or excessive inhomogeneties in structures 
such as radome structures, and other applications. 

BRIEF SUMMARY 

Therefore, it is a primary object, feature, or advantage of 
the present invention to improve over the state of the art. 

It is a further object, feature, or advantage of the present 
invention to provide a sensor and related systems and meth- 
ods which may be used in quantitative characterization of 
material properties of multi-layered planar dielectric struc- 
tures. 

Yet a further object, feature, or advantage of the present 
invention is to provide a sensor and related systems and 
methods for use in detecting water or excessive inhomogene- 
ities caused by repairs in modem radome structures. 

A still further object, feature, or advantage of the present 
invention is to provide a sensor and related systems and 
methods which may be used in quantitative characterization 
of material properties of multi-layered cylindrical dielectric 
structures. 

Another object, feature, or advantage of the present inven- 
tion is to provide a sensor and related systems and methods 
which are appropriate for use in handheld devices. 

Yet another object, feature, or advantage of the present 
invention is to provide a rotationally-invariant capacitive 
probe. 

A still further object, feature, or advantage of the present 
invention is to provide for capacitive probes that allow 
removal of parasitic capacitances. 

The present invention is not to be limited to or by these 
objects, features, and advantages. It is to be further under- 
stood that no single embodiment of the present invention need 
exhibit all of these objects, features, or advantages. 

According to one aspect of the present invention, a con- 
centric coplanar capacitive sensor is provided. The sensor 
includes a charged central disc fomiing a first electrode and 
an outer annular ring coplanar with and outer to the charged 
central disc, the outer annular ring forming a second elec- 
trode. There is a gap between the charged central disc and the 
outer annular ring. The first electrode and the second elec- 
trode are attached to an insulative film. 

According to another aspect of the present invention, a 
capacitive nondestructive evaluation system for evaluating a 
dielectric test piece is provided. The system includes a con- 
centric coplanar capacitive sensor having (a) a charged cen- 
tral disc forming a first electrode, (b) an outer annular ring 
coplanar with and outer to the charged central disc, the outer 
annular ring forming a second electrode, and (c) a gap 
between the charged central disc and the outer annular ring. 
The system also includes a capacitance measuring circuit 
electrically connected to the concentric coplanar capacitive 
sensor for measuring transcapacitance between the first elec- 
trode and the second electrode for use in evaluating the dielec- 
tric test piece. 

According to another aspect of the present invention, a 
method of non-destructive evaluation is provided. The 
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method includes providing a concentric coplanar capacitive 
sensor, attaching the concentric coplanar capacitive sensor to 
a dielectric test piece, applying an input signal across the 
concentric coplanar capacitive sensor to produce an output 
5 signal, determining transcapacitance between the first elec- 
trode and the second electrode based on the output signal, and 
using the transcapacitance in a model that accounts for the 
dielectric test piece to determine inversely the properties of 

the dielectric test piece, 
to 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a concentric coplanar capacitive sensor. 
The radius of the central disc and the width of the outer ring 
are denoted s and t, respectively. The gap between them is g, 
and D is the sensor diameter 

FIG. 2 illustrates a point charge on top of a four-layer 
dielectric. 

, 0 FIG. 3 illustrates a point charge on top of a layered half- 
space dielectric. 

FIG. 4 illustrates a concentric sensor divided into N circu- 
lar filaments, each with a constant surface charge density that 
is constant with respect to variation in p 
25 FIG. 5 illustrates a calculated surface charge distribution 
for the sensor shown in FIG. 1 in contact with a half-space 
dielectric. Sensor configuration: s=t=10 mm, g=l mm, 
V„„, ar = 1 V, and V OMer =0 V. The test-piece has relative dielec- 
tric constant E,.=8. 

30 FIG. 6 is a graph illustrating a normalized sensor tran- 
scapacitance versus changing sensor disc radius s and elec- 
trode gap g. The sensor outer radius D/2=s+g+t is fixed. 

FIG. 7 illustrates how a calculated sensor output signal IC r l 
changes as a function of E rl =E r3 and the core-layer relative 
permittivity E r2 . IC r l is normalized by its own maximum 
value for this calculation, which is 4.66 pF. Sensor configu- 
ration is as for FIG. 5. 

FIG. 8 illustrates measured and calculated \C T \ for various 
40 sensor configurations (see FIG. 1) in contact with a glass plate 
with E,.=5.62 and thickness 3.02 mm. 

FIG. 9 illustrates measured and calculated IC r l for various 
sensor configurations (see FIG. 1) in contact with an acrylic 
plate, E r =2.85 and thickness 2.39 mm, on top of a glass plate 
45 with parameters as for FIG. 8. 

FIG. 10 illustrates measured and calculated \C T \ for vari- 
ous sensor configurations (see FIG. 1) in contact with a three- 
layer acrylic-glass-acrylic structure. Layer parameters are as 
for FIGS. 8 and 9. 

50 FIG. 11 illustrates a sensor on top of a 1 cc water- injected 
glassfiber-honeycomb-glassfiber sandwich panel. The sub- 
figure is a photograph of the sandwich panel whose properties 
are given in Table 1 . 

FIG. 12 illustrates measured IC r l for 1 cc of water injected 
55 into the glassfiber-honeycomb-glassfiber sandwich panel, 
FIG. 11. Sensor configuration is as for FIG. 5. 

FIG. 13 is a block diagram illustrating one example of a 
capacitive NDE system which uses the sensor. 

FIGS. 14A and 14B provide illustrations of concentric 
60 capacitive electrodes on top of a multi-layer dielectrics: (a) 
sensor configuration and test-piece structure used in the 
numerical modeling; (b) assembled hand-held probe based on 
the modeling in (a). 

FIGS. 15A and 15B provide photograph of the assembled 
65 probe: FIG. 15A is an experiment setup used in probe lift-off 
measurements; FIG. 15B is a concentric capacitive sensor 
fabricated by photolithography. 
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FIG. 16 is a block diagrams of the assembled hand-held 
probes illustrating equipment used in capacitance measure- 
ments and equipment used in probe calibration. 

FIGS. 17A-17C penetration depth of concentric capacitive 
sensors. FIG. 17A is a cross-section view of a concentric 5 
capacitive sensor in surface contact with a one-layer dielec- 
tric in free space. FIG. 17B is a calculated sensor output 
capacitance as a function of test -piece permittivity and thick- 
ness. FIG. 17C is a calculated difference between the capaci- 
tance in FIG. 17B and that of a similar but infinitely thick 
test-piece. 

FIGS. 18A-18D illustrate measured capacitance of hand- 
held probes as a function of test-piece thickness. FIG. 18A 
illustrates stepped Delrin® slab £,.=3.82; FIG. 18B illustrates 
stepped HDPE slab G r =2.65. FIGS. 18C and 18D illustrate 
difference calculated for Delrin® and HDPE respectively, 15 
assuming an 18-mm-thick test-piece to be an approximate 
half-space 

FIGS. 19A-19D measured and calculated differences in 
capacitance of hand-held probes as a function of liftoff. FIG. 
19A is for sensor A. FIG. 19B is for sensor B. FIG. 19C and 20 
FIG. 19D illustrate permittivity determined for PMMA and 
glass, respectively. 

FIG. 20 is a photograph of the sandwich panel with param- 
eters listed in Table 4. 

FIGS. 21A-21D illustrate capacitance measured as hand- 
held probes scan over glassfiber-honeycomb-glassfiber sand- 
wich panels containing injected dielectric contrast agents 
water and olive oil. FIG. 21A is sensor A and injected water. 
FIG. 21B is sensor A and injected olive oil. FIG. 21C is sensor 
B and injected water. FIG. 21D is sensor B and injected olive 
oil. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

1 . Introduction 35 


shown having a charged central disc 14. There is an outer 
annular ring 12 which is coplanar with and outer to the 
charged central disc 14. There is a gap 1 6 between the charged 
central disc 14 and the outer annular ring 12. 

Electrostatic Green’s functions due to a point charge over 
different test-piece structures are derived first. These Green’ s 
functions are then utilized in later MoM calculations of the 
sensor transcapacitance C T . Because of the cylindrical sym- 
metry of the designed sensor, the electrostatic Green’s func- 
tions are derived in cylindrical coordinates through the Han- 
kel transform method. Additionally, the test-pieces in our 
theoretical analyses are assumed to be infinite in the horizon- 
tal directions and the sensor electrodes are assumed to be 
infinitesimally thin. 

Assume there is a point charge placed at the origin in free 
space. The resulting electrostatic potential 'P, related to the 
electric field E=-V*P, satisfies the Laplace equation and can 
be expressed in cylindrical coordinates as 


< a 2 id a 2 \ (i) 

lv + P^ + 4?r p ’ i)=a ''*° 

where T^p, z) is independent of azimuthal angle j). Next, the 
Hankel transform ;f(K) of zero-order of a function ;f(p) is 
given by 

(2) 

where J 0 (z) is the Bessel function of the first kind and the 
inverse transform is of the same form. Apply the zero-order 
Hankel transform to equation (1 ), making use of the following 
identity 


d 2 l d 
dp 2 p pdp 


\f(p)\j 0 (Kp)pdp = -K 2 f(K), 


( 3 ) 


The present invention provides for a concentric coplanar 
capacitive sensor which may be used for detecting water or 
excessive inhomogeneities caused by repairs in modern 
radome structures. The proposed sensor, having the advan- 
tage of rotational symmetry, consists of a charged central disc 
and a coplanar outer annular ring that exhibit a measurable 
transcapacitance C T . The output signal depends on the mate- 
rial and structural properties of the test -piece with which the 
sensor is in surface contact. An electrostatic Green’s function 
for a three-layered dielectric structure in free space is derived 
in cylindrical coordinates through the Hankel transform 
method. This derived Green’s function may then be simpli- 
fied, providing results for many other cases such as a half- 
space dielectric, a layered half-space dielectric, and one- and 
two-layered dielectrics in free space. Numerical implemen- 
tations based on these Green’s functions are described, in 
which the surface charge distribution on the sensor electrodes 
is calculated through the method of moments (MoM). From 
the surface charge, C r is calculated. To verify the validity of 
the numerical calculation, benchmark experiments are con- 
ducted for one-, two-, and three-layer dielectric test-pieces in 
free space, respectively. Very good agreement is observed 
between the calculated and measured transcapacitance. Fur- 
thermore, water ingression measurements in a sandwich 
structure are carried out and demonstrate the feasibility of 
using the capacitive sensor to detect water intrusion and inho- 
mogeneities in radome structures. 


where ;f(p) is assumed to be such that the terms pJ 0 (Kp)3;f 
(p)/3p and pf (p)3J 0 (Kp)/3p vanish at both limits. The spatial 
domain Laplace equation (1) is then transformed into a one- 
dimensional Helmholtz equation in the transformed domain: 


( & , 1 - ( 4 ) 

where for k the root with positive real part is taken. Here, the 
Green’s functions are first derived in the transformed domain 
and then transformed back to the spatial domain through the 
inverse Hankel transform. 

The present invention further contemplates that sensor sur- 
face charge density may be computed in alternative ways. For 
example, the spectral domain Green’s function may be used 
to derive the integral equation for the sensor surface charge 
density in the spectral domain, using Parseval’s theorem. 
Then the integral equation may be discretized to form matrix 
equations using the MoM. The spatial domain approach is 
more computationally efficient for both one- and three-lay- 
ered structures in free space, while the Green’s function deri- 
vation and numerical implementation for the spectral domain 
approach are more straightforward. Additional details regard- 
ing the alternative approach are described in [24]. 


40 


45 


50 


55 


60 


2. Green’s Functions for Multilayered Dielectrics 


2.1 Point Chaige on Top of a Four-Layer Dielectric 


The configuration of the proposed sensor is shown in FIG. 
1. In FIG. 1, a concentric coplanar capacitive sensor 10 is 


65 


One potential application of the capacitive sensor is the 
dielectric property characterization of three-layer modern air- 
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craft radome structures, using the knowledge of sensor geom- 
etry and the output transcapacitance C T . In order to set up the 
governing equations in the MoM calculations for the in- 
contact characterization of layered dielectric structures, the 
potential due to a point charge in the plane z=0 is derived. 5 
Without loss of generality, a four-layer half-space dielectric 
configuration shown in FIG. 2 is used in the following theo- 
retical derivation. One can easily obtain the solution for the 
three-layer radome structure by replacing <E 4 by (E 0 , the per- 
mittivity of free space. to 

In FIG. 2, a point charge is placed on top of a four-layer 
half-space dielectric. The electrostatic potential *P satisfies 
the Laplace equation in each homogeneous medium. After 
applying the zero-order Hankel transform mentioned above, 
the resulting one-dimensional Helmholtz equations in the 15 
transformed domain can be expressed as 


Z > 0 

(J-* 2 K'Z) = 0 

- h - : < — A/_i 


(5) 


20 


( 6 ) 


25 


where i=l , 2, 3, 4, and ho=0 while h 4 -»- oo. The subscripts 0, 

1 , . . . , 4 denote the free space above the dielectric and each 
homogeneous layer of the dielectric, respectively. From equa- 
tions (5) and (6), general solutions for the potentials in each 30 
region can be expressed as 

'V t (K^)=A i (K)e~ KZ +B l (K)e KI , -hjSz<-h,_i, (7) 

where B o (k)=A 4 (k)= 0 due to the fact that the potential at 
infinity vanishes. 35 

The interface conditions on the electric fields are 


8 

Therefore, the surface charge density in the plane z=0 is 
p s =S(p)/2jtp, with its Hankel transform being Vat. Applying 
the Hankel transform to the boundary condition equation (8), 
one can easily get the result shown in equation (11). 

Substitute equation (7) into equations (10) to (13) to 
express the coefficient Aj(k) as 


Ai (*) = 


(15) 


x L«5e -2 ** 3 + ye-™'!* pe- lKh ' +/Jy < 5e~ 2 «*i + ' 1 3-'>2)j * 

Z71K{Sq +£ i ) 


Z 

71=0 


HD" 


(/? + are -2 * 7 ! )” 
(i+a/3e- 2l,T if +l 


(y + de -2 "^)" _. 

(1 + y6e~ 2l<T 3) n+1 £ 


where a=(E 1 -E 0 )/(e 1 +e 0 ), p=(e 2 -e 1 )/(e 2 +e 1 ), Y=(G 3 - 
E 2 )/(<E 3 +<E 2 ), 8=(eE 4 -E 3 )/(E 4 +e 3 ). Besides, T^h,, T 2 =h 2 - 
h 15 and T 3 =h 3 -h 2 . In order to get the spatial domain solution, 
equation (15) can be expanded into the form of series sum- 
mations, which facilitates application of the inverse Hankel 
transform. For those terms inside the summation of A x (k), we 
have 




r\{n-r)\ 


/?-Ve- 2 " T i , 


i = f 

( 1 + a/3e~ 2 * T i ) n+l «!^! 


(16) 


(17) 


and similarly for terms (y+be- 2 ' cT, y and (l+y8e- 2Kr3 )- (,,+1) . 
Combining equations (16) and (17) gives 


fxfEo-E^O, £-(p 0 -D,)= Ps 
fx(£,.-£ j+1 )=0,f-(A-A + i)=0 


( 8 ) 


(9) 40 


{ p +ae -^f (18) 


(1 + crpe-^i )" 


r\(n - r)!.s! 


where 1 = 1 , 2, 3, and p is the free surface charge density in the , . , , . . , . „ 

, nti- tt i i ■ 7 Adopting m=r+s, equation (18) is then written in the follow- 

plane z=0. Applying the Hankel transform to the interface . 2 ^ n ' 

conditions for E and D, the corresponding boundary condi- ln ® 0rm 

tions for the potentials in the transformed domain are 


expressed: 


45 




( 1 + a fie : 

%(k, o yM%(K, 0 ), 

(10) 

where 

d'Pofo 0) (at, 0) 1 

£ ° dz 1 CL dz 2tt ’ 

(11) 50 

K mn (a, P) 

%(k, -hi)=y i+l (K, -hi). 

(12) 


d^i(K, -hi) d^i+i(K, -hi) 

(13) 

Similarly 


. ,71+1 “ Zj 
> 771=0 


(m + n — r ) ! 
r\(m -r)\{n-r) 


oTp m+ ‘ 


(19) 


( 20 ) 


dz 


dz 


where i=l , 2, 3. A little more explanation is made here about 
the Vat term on the right-hand side of equation (11). In cylin- 
drical coordinates, the Dirac delta-function can be expressed 60 

for points on the z axis as „. „ , . ...... 

Finally, the senes summation fonn lor A^k) in the trans- 
formed domain is written as 


(y + <Se _2 "3)" 
(1 +y6c- 2 * 7 '3)" 4 


= Y J K ln {S,y)e- mT 3. 


( 21 ) 


65 


d(r -/)= ^— 6(p)6(z ~ Z ') 

Inp 


( 14 ) 


4iW = 


(22) 
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-continued 

[<5e- M 3 +ye -2«*2 + y Syd>- 2 "<*l + *3-*2)] X 

2 jIK(£q +£l) 

Z Z Z j8)*fa(<5. ric -2 '" 7 ' 1 e-^e- 2 ^ , 

w=0 /=0 »i=0 

and it is found from the boundary conditions that 


10 

throughout the entire domain can be derived from the above 
equations but only ^(p.O) is needed here for later MoM 
calculations because the sensor is in contact with the test- 
piece surface. By substituting E 0 for E 4 in the above relations, 
5 the potential due to a point charge on top of a three-layer 
dielectric in free space can be retrieved. Numerical results 
based on this potential are compared with corresponding 
experimental results in Section 4. 


/toW 


1 

2tCK(£q + £\ ) 


+ (1 + o)A\{k) 


10 

(23) 


2.2 Point Charge on Top of a Two-Layer Dielectric 
in Free Space 


Applying the inverse Hankel transform to equation (23), the 15 
potential in the z=0 plane due to a point charge at the origin is 
expressed as 


1 (24) 20 

%(p,0)=— -+(1+0) Al(K)Jo(Kp)KdK. 

2tt(£o +£i)p Jo 


The integral in equation (24) can be evaluated by applying the 
following Hankel transform pair to each of its power series 
terms, given in equation (22), 


25 


r°°e " 1 

I Jn(Kp)KdK- 

J 0 K 


(25) 


Vp 2 + z 2 

Equation (24) is finally expressed in real-space form as 


30 


1 1 + a (26) 

Vo(AO)=— — - — — - — -( Gl +G2 + G 3 +G 4 ), 


35 


27T(£ 0 + £i )p 2 n(£ 0 + £i ) 


where 


°.-ZZZ (-1 fK m [a,P) 

n = 0 1=0 m = 0 

K,„(S,y)x. 


(27) 


ojp 2 + PtmT, + nT= + IT} + A 3 )] 2 


Gi = ZEE (-1 TK mn (a,P) 

n = 0 1=0 m = 0 

Kin(S, y)x 


(28) 


40 


45 


°“EEE (-1 r KUa.fi) 

n = 0 /=0 m = 0 


ojp 2 + [2(mT l +nT 2 + 1T } + /i 2 )] 2 5Q 

(29) 


A' fc (Ay)x 




Vp 2 


+ [2(mTi +nT 2 +lT 3 + hi )] 2 


G 4 = EEE (-lfA m „(o./J)A' to ((5.y)x 

n =0 1=0 m =0 

PyS 


(30) 


'J p 2 + + //'A + /73 + A, + /+ — A 2)] 2 


Equations (26) to (30) together give the surface potential 65 
'Pq(p.O) due to a point charge at the surface of a four-layer 
half-space dielectric in the spatial domain. The potential 


The surface potential for the case of a point charge on top 
of a two-layer dielectric can be simplified from equation (26) 
by assuming that Ez l =€z 1 and E 4 =€E 0 . We are interested in this 
case for the purpose of benchmark testing described in Sec- 
tion 4. As a result, (3 becomes zero and G 3 =G 4 =0. Ontheother 
hand, K m „(a,(3) lias a non-zero value, K mn (a)=a’\ only when 
m=n=r. This is because when m*n, the term m+n-2r is con- 
stantly greater than zero and thus P" !+ ” _2 ' =0. Hence, the cor- 
responding potential is simplified as 


%(P. 0) = 


1 

2n(so +£i)p 


1 + a 

2ir(£o +£1) 


(G 1 + G 2 ), 


where 


(31) 


(32) 


yy (-ifa , K.g,r)x- =. 

^0 7^0 Vp 2 + {2[(n + l)7' 1 +(m + l)r 2 ]} 2 




(33) 


G 2 = > > (-l)VA m „(<5, y)x- 

n =o ,„= 0 a/p 2 + 12[(« + 1)7-1 +mr 2 ]} 2 


Here, Tj and T 2 represent the thickness of the top and bottom 
homogeneous layers, respectively. 


2.3 Point Charge on Top of a Two-Layer Half-Space 
Dielectric 


The above derived potential due to a point charge over the 
surface of a four-layer half-space dielectric can also be 
reduced to the case of a point charge on top of a coated 
half-space dielectric. This case can be furthermore reduced to 
the solutions of a point charge on top of a one-layer dielectric 
slab in free space and a point charge on top of a homogeneous 
half-space dielectric. These simplified results are identical to 
those presented in [21] and [22]. In addition, calculation 
results based on the potential due to a point charge on top of 
a one-layer dielectric in free space are used in the benchmark 
comparison in Section 4. 

Assuming that E 1 =E 2 =E 3 *E 0 , the structure in FIG. 2 is 
simplified into the case of a half-space dielectric with a single 
surface layer as shown in FIG. 3. The top layer has dielectric 
constant E x and thickness h. The bottom layer is the half- 
space dielectric with dielectric constant E 2 . In this case, 
P=y=0. K m „(a,[3) only has non-zero value when m=n=r and 
K„,„(a)=a". Similarly, K ? „(8,y) only has non-zero value when 
l=n=t and kj„( 8)=8". Equation (26) is simplified to 
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%(p, 0) = 


(34) 


1 

2x(so +£l) 


x \ - - (1 + a) 


z 


(- 1 )” 


(af(6r 


+J p* +[2(n+-l)h] 2 


5 


'P(p, <p,0\p',tp',0) = 


2^(fio 


1 1 

) +«i) X \|r-r'| 


-(!+<*) 


Gi(|p-p'|) + G,(|p-p'|)+j) 

G 3 (lp-p / |) + G 4 (|p-p'|) J| 


when? 


( 37 ) 


\r — r f \ = Vp 2 + p 2 — 3p//cos(7' — </>'), 


( 38 ) 


where a=(G 1 -G o )/(G 1 +0 o ) and 8=(E 2 -8 1 )/(E 2 +E 1 ). To 
compare the derived result with that in the literature, rewrite 10 
equation (34) in terms of coefficients a=(G 1 -e o )/(G 1 +e o ) 
and (3=(S 1 -E 2 )/(E 2 +Ei), which gives 

15 

%(p, 0) = (35) 


Gi(|p-p'|) = 


ZEE 


(-If KU«, P) 


( 39 ) 


k,„(S, r) x 

Vk-P I + [ 2 (m 7 i +11T2 +IT) + /13)] 2 


1 


2 n(s 0 +£i) 


xl -+(! + «)> (apf 


z> 


Vp 2 + [2(n + l)/l] 2 


20 


Equation (35) is identical with the result presented in [21], 
where the Green’s function is derived using a double Fourier 
transform in Cartesian coordinates. 

A special case is that in which the half-space dielectric is 
replaced by free space and the test-piece in contact with the 
sensor is then a homogeneous plate. The corresponding 
potential is expressed in equation (36) by replacing €E 2 with 
E 0 in equation (35): 


and G,(lp-p’l), G 3 (lp— p’l), and G 4 (lp-p'l) can be modified 
similarly. For other test -piece configurations, the appropriate 
Green’s function should be used. Moreover, the potential at 
such an observation point due to points on a charged sensor 
shown in FIG. 1 can be derived by integrating equation (37) 
over the sensor electrode surface: 

t(p. 0|p', 0) = 1 — -x f K(p,0\p\0)cr(p’)p'dp\ 

2?T(£o+«l) Jdisc+ring 


where a(p') is the sensor surface charge density and 


%(P, 0) = 


(36) 


X i - + (1 + a) 

2x(£ 0 +£i) |p 


00 ^ \ 35 r 4 

V . A-(p.0|p'.0) = f -J—-(l +a )YG,(lp-p'l) 


W- 


(41) 


Equation (36) can be simplified further by choosing 
G^Eo- The series summation terms in equation (36) all van- 
ish because a=0 in this case. This simplified result is identical 
to that presented in [22], in which the result is derived in the 
spatial domain directly. 

3. Numerical Implementation 
3.1 Calculation Method 

In order to calculate the sensor transcapacitance, Cj, the 
method of moments (MoM) [23] is utilized in the numerical 
calculations. In the following calculation examples, all the 
sensors share the configuration shown in FIG. 1, where the 
central disc is charged to the potential Vj=l V and potential of 
the outer ring is kept at V 2 =0 V. 

The electrostatic potentials due to a point source, derived 
above, serve as the Green’s functions in the MoM simula- 
tions. As shown in FIG. 4 , the concentric electrodes of the 
sensor are divided into N circular filaments each with width A 
and a surface charge density that is constant with respect to 
variation in p. For the test-piece structure shown in FIG. 2 , the 
potential at a given observation point (p, t[), 0) due to a source 
point (p’, 4>’, 0) can be expressed as follows, by slightly modi- 
fying equations (26) to (30): 


One thing to notice is that because of the cylindrical sym- 
metry of the sensor structure, the resulting potential in space 
is independent of the azimuthal angle (|>. Therefore, the prob- 
lem of calculating the sensor surface charge distribution, 
which is determined by the potential distribution, is reduced 
to the p -direction only. For observation points on the sensor 
electrodes, the boundary conditions for the potential can be 
expressed as 

%(p, z = 0)= 1 x f K{p, 0 1 p\ 0)o-(p')p’dp' = V m , <4-) 

2tt(£o + £l ) Jdisc+ring 

where points on the central disc are denoted by m=l while 
those on the outer ring are denoted by m=2. In order to solve 
for the sensor surface charge distribution cr(p') using MoM 
calculations, the following expansion is used: 

n ( 43 ) 

°V) = Yj °"AV) 

j = 1 

where tr (p') is the basis function and Oj is the unknown 
coefficient. Here, we choose b 7 (p 1 ) as the following function 
for filaments on the inner disc 


40 
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r i 

bj ip') = | -Js 2 - (p') 2 
[ 0, 


(y-l)A<p' <;A 


elsewhere. 


(44) 


5 


where s is the radius of the inner disc. For the filaments on the 
outer annular ring, b .(p 1 ) is chosen as 

to 

bj(p') = < 45 ) 

l 1 1 

, X— = O’- 1)A <p' < jA 

I Vo + ?) 2 -(P'l 2 y(D/2) 2 -(p') 2 

( 0, elsewhere, ^ 


where g is the gap between the two sensor electrodes and D is 
the diameter of the sensor. This fonn of basis function has the 
advantage of modeling the edge effect of the charge distribu- 2 0 
tion discussed later. To resolve the N unknown Oj coefficients, 
it is then required that the boundary conditions for V m in 
equation (42) are satisfied for each circular filament on the 
sensor surface. To evaluate equation (42) in N different fila- 
ments, weighting (or testing) functions w,.(p) are needed. 
Here, we choose the weighting and basis functions to be the 
same, known as Galerkin’s method. For filaments on the inner 
disc, 


W;(p) = 


1 

V^-tp') 2 

0, 


(i- 1)A <p' < /A 
elsewhere, 


(46) 


35 


while the weighting function for filaments on the outer annu- 
lar ring 


40 

(1 1 , (47) 

— x — (i - 1)A < p < /A 

wi(p) = \ y (s + g ) 2 - (pO 2 ^ (D/2) 2 - ip') 2 

{ 0, elsewhere, 

45 

where i=l, 2, . . . , N. Discretizing the integral equation using 
weighting functions in each of the N filaments, equation (42) 
turns into the following matrix equation: 


Gu 

G12 

■ Guv 


<r 1 ‘ 

G21 

G22 

• G 2 a/ 

X 

2 

. GWi 

G/V22 ■ 

• &nn . 


o' N . 


where 


50 

(48) 


55 


po-m po-DA 

Gj= I w,- (p)x I A'(p. 0|p', 0)6 y (p')p' dp’ pdp. 
JiA LU jA 


(49) 


For the V matrix, if the element is located on the central 
electrode, its value is V 1 =l V; while the values for those 
elements located on the outer ring are V 2 =0. 

From equation (48), the sensor surface charge distribution 
CT(p') can be calculated. Once cr(p') is known, one can inte- 65 
grate over the electrode surfaces and find the total chaige on 
both imier and outer electrodes. The sensor output signal, 


14 

which is the transcapacitance C T between those two elec- 
trodes, can be ultimately calculated through 


c _ Qouter I (50) 

Vinner Wouter =(> 

where Q ollter is the total charge on the outer electrode, while 
V in „ er and V oute; . respectively represent the voltage on the 
imier and outer electrodes. Choosing this convention leads to 
C r <0, whereas ICjl is compared with experiment in the fol- 
lowing. 


3.2 Example Calculations 

FIG. 5 shows an example of the sensor surface charge 
distribution, where the sensor is placed above a half-space 
dielectric with relative permittivity E,.=8. The sensor configu- 
ration is s=t=10 mm and g=l mm Due to the edge effect, the 
surface charge density at the edge of the imier charged elec- 
trode is singular. This positive charge distribution results in a 
negative surface charge distribution on the outer electrode. 
The surface charge density on the imier edge of the outer 
electrode tends to infinity much faster than that on the outer 
edge, because of its smaller radius and stronger interaction 
with the imier electrode. It is worth mentioning that when one 
applies a different combination of potentials on the inner and 
outer electrodes, the sensor surface charge distribution 
changes correspondingly. However, the sensor transcapaci- 
tance ICjl, which is the intrinsic property of the sensor and 
only determined by its own structure, is unchanged. The 
sensor transcapacitance, IC r l=5.398 pF for this case, is cal- 
culated through equation (50). 

Numerical calculations based on the same test-piece have 
been carried out to investigate the optimal sensor configura- 
tion giving the maximum output signal IC r l. The sensor out- 
put signal as a function of s and g is plotted in FIG. 6. In the 
calculation, the sensor outer radius D/2=s+g+t is fixed and all 
the curves in FIG. 6 are normalized with respect to their own 
maximum values. As can be seen from the figure, for any 
given g, the sensor output signal increases to a maximum 
value and then decreases as s increases. This is because as s 
increases, the width of the outer electrode t decreases, result- 
ing in stronger edge effects on its surface charge distribution. 
These stronger edge effects result in more charges accumu- 
lated on the outer ring, and therefore the sensor output signal 
is increased according to equation (48). In this regime, the 
surface charge density is the dominant factor determining the 
total surface charge Q outer . However, as s increases and passes 
a certain value, the sensor output signal starts to decrease. 
This is due to the fact that the diminishing surface area of the 
outer electrode becomes dominant in determining the total 
surface chaige Q outer . As a result, we observe an optimal 
sensor configuration for a given g that gives the maximum 
ICyl . It is also verified in our calculations that the shape of all 
the curves in FIG. 6 do not depend on the actual size of the 
sensor and the applied electrode voltage, but only on the 
relative values of s, g, and D. Similarly, as g increases, the 
interaction between the imier and outer electrodes is 
decreased, and the surface charge density at their neighboring 
edges diminishes accordingly. Because of the decreased edge 
effect and surface charge density, the outer electrode needs 
more surface area to achieve its maximum Q outer , which is 
directly proportional to ICjl . Hiis is why as g/D increases, the 
s/D value that yields the maximum IC r l decreases in FIG. 6. 
As one can imagine, the absolute magnitude of ICy also 
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becomes smaller for larger g and fixed s and D values, due to 
the same reasoning mentioned above. Consequently, in order 
to achieve the maximum IC r l, it is desirable to maintain high 
s/D and low g/D ratios. Nevertheless, it is worth mentioning 
that the sensitive area of the sensor closely corresponds to the 5 
location of the gap between its two electrodes, and there will 
be an insensitive zone at the center of those sensors with 
relatively large s values. 

Another example, addressing sensor sensitivity to changes 
in core permittivity of a three-layer structure, is presented 
here. We are interested in this problem because one potential 
application of the sensor is detection of ingressed water or 
inhomogeneities in the core of an aircraft radome structure, 
which is typically a three-layer sandwich structure. In the 
numerical calculation, the infinite series summations in equa- 
tions (27) to (30) are truncated to N=10 terms each. The 15 
difference between N=10 and N=100 terms is only 0.008% 
while the latter is extremely time-consuming. The sensor 
configuration is s=t=10 mm and g=0.5 mm. The test-piece is 
shown in FIG. 2, where T 1 =T 3 =24 mm, T 2 =3 mm, and 
medium 4 is replaced by free space. The relative permittivity 20 
of the top and bottom layers, Eq and E 3 , is chosen to be the 
same. These parameters are also adopted in later benchmark 
experiments described in Section 4. FIG. 7 shows how the 
normalized sensor output signal IC r l changes as a function of 
Ej=E 3 and of the core relative permittivity E 2 . In FIG. 7 , 25 
E rl =E r2 =E r3 =l gives the limiting case of the sensor in free 
space; E rl =E r2 =E r3 ^l gives the case of the sensor on top of 
a one-layer test-piece in free space; and E rl =E,. 3 = 1 *E 2 gives 
the case of lift-off measurement of a one-layer test-piece in 
free space. It is seen from FIG. 7 that the slope of the curve 30 
representing the normalized IC r l as a result of changing 
E rl =E. r3 when E . r2 = 1 0 is much greater than that obtained as a 
result of changing E,. 2 when E rl =E,. 3 =10 as expected due to 
the shielding effect of the top layer. In addition, high E rl =E,, 3 
values give less sensitivity to e r2 changes. This can be made 35 
more explicit by defining the percentage difference in the 
sensor output signal as follows: 


% difference = P = — x 100%. 

\Ct\^ 


(51) 40 


When E rl =E r3 =3, E,. 2 =2, and AE,. 2 =1 , for example, then P is 
3.66%. However, for the same E,, 2 and AE,. 2 , when 45 
E,. 1 =E,. 3 =10, P is only 2.99%. This percentage change in IC r l 
is expected to be even smaller when E rl becomes larger, 
which is reasonable because higher density electric fields are 
confined in the high E rl material. To improve sensor sensi- 
tivity to the permittivity change in the core-layer then, one can 50 
increase the gap g between the electrodes to some extent. For 
example, when g=l mm rather than 0.5 nun as in the calcu- 
lations of FIG. 7, and keeping all the other parameters the 
same, P is 3.62% when E rl =E r3 =3 and 4.42% when 
E rl =E r3 =10. However, the magnitude of the sensor output 55 
signal is decreased as g increases. Therefore, a trade-off 
between high sensor sensitivity and strong output signal is 
needed when determining the optimal sensor configuration 
for measurements detecting permittivity change in the core 
layer. 60 


4. Experiments 


4.1 Benchmark Experiments 

In order to verify the validity of the theory developed 
above, benchmark experiments were carried out for one-, 
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two, and three-layer dielectric test-pieces in free space, 
respectively. An Agilent E4980A precision LCR meter (20 
Hz to 2 MHz) was utilized for the capacitance measurements . 
The operating frequency of the LCR meter was set to be 1 
MHz. This particular frequency ensured that the measure- 
ment error of the LCR meter was less than 0.3% for a 1 pF 
capacitance, while at the same time giving a good approxi- 
mation for the electrostatic case in the numerical model. A 
Novocontrol Alpha Dielectric Spectrometer was used to inde- 
pendently measure the dielectric constants of the samples 
used in the benchmark experiments . In the Novocontrol mea- 
surements, two 40-nun-diameter electrodes were used and 
the edge effect compensation was turned on, due to the fact 
that the thicknesses of the test-pieces were relatively large 
compared to the test fixture’s electrode diameter. In addition, 
the test-piece thicknesses were measured by a digital thick- 
ness indicator with accuracy ±1 pm. These independently- 
measured test-piece thickness and dielectric constant values 
were used as the inputs of the calculation model. 

Seven copper sensors of the configuration shown in FIG. 1 
were fabricated by photolithography. Four sensors have 
g=0.5 nun and three have g=l nun, with different s=t values. 
These sensors were deposited onto a 25-pm-thick Kapton® 
film to support the copper. By comparing the calculation 
result of a capacitive sensor (s=t=10 mm and g=0.5 mm) on 
top of a half-space dielectric (E,.=8) and that of the same 
sensor on top of a 25-pm-thick Kapton® film over the same 
half-space, it was estimated that the presence of the Kapton® 
film influences the measurement signal by less than 0.5%. For 
each of the following benchmark measurements, the test- 
piece was supported by three acrylic stands 50 cm above a 
wood-top working table to approximate the infinite test -piece 
in free space assumption in the calculation model. Tape was 
used to attach each sensor tightly against the test-piece to 
ensure minimum air gap between the sensor and the test- 
piece, due to the fact that the presence of an air gap can affect 
measurement results significantly. The tape was attached on 
the edges of the Kapton® film, far away from the sensor outer 
electrode. IC r l was measured by placing the probe of the 
Agilent probe test fixture 1 6095A across the two sensor elec- 
trodes. This probe test fixture was comiected to the LCR 
meter and the capacitance values were read from the LCR 
screen. 

To verify the results for the case of the capacitive sensor on 
top of a one-layer dielectric test-piece in free space, a glass 
plate with dimensions 305x305 nun 2 and thickness 3.02±0.01 
nun was used. The test -piece dielectric constant was indepen- 
dently measured as 5.62±0.05. FIG. 8 gives the comparison 
between the numerical and experimental results. Experimen- 
tal data show excellent agreement with the numerical results, 
to within 4%. Ten measurements were made for each sensor 
and the results were averaged. The maximum standard devia- 
tion in the measurements was found to be 2% . As can be seen, 
ICyl increases as s increases and decreases as g increases. 
Meanwhile, sensors with smaller s values show relatively 
greater standard deviation in the measured data. This is rea- 
sonable because when the scale of the sensor becomes 
smaller, the output capacitance is consequently smaller, and 
the noise from the surroundings in the measurement environ- 
ment can have a relatively greater impact on the measurement 
results. 

The case of the capacitive sensor on top of a two-layer 
dielectric test-piece in free space was verified by placing a 
305x305 nun 2 acrylic plate with thickness 2.39±0.02 nun on 
top of the glass plate mentioned above. The independently 
measured acrylic dielectric constant was 2.85±0.05 in this 
case. Plastic clamps were used to make sure there was as little 
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air gap as possible between these two plates. FIG. 9 gives the 
comparison between the numerical and experimental results. 
Again, very good agreement between experimental and theo- 
retical results is observed. The maximum difference between 
the theory and experiment is less than 3% and the maximum 5 
standard deviation is 1% in these measurements. Similarly, 
FIG. 10 shows the comparison results for the case of the 
capacitive sensor on top of a three-layered acrylic-glass- 
acrylic structure. The top and bottom acrylic plates share the 
same parameters and the glass plate sandwiched in the middle to 
is the same as that used previously. It is seen from FIG. 10 
that, even for this more complex test-piece, very good agree- 
ment between theoretical predictions and experimental 
results is obtained. In this case, the maximum difference 
between the theory and experiment is 3% and the maximum 15 
standard deviation is 1%. 

In conclusion, benchmark experiments show very good 
agreement with theoretical predictions. The output signal for 
the three-layer acrylic-glass-acrylic structure is slightly 
greater than that of the two-layer acrylic-glass structure but 20 
smaller than that of the one-layer glass plate. Because glass 
has a higher permittivity than acrylic, the sensor output signal 
of the one-layer glass plate is greater than that of the two-layer 
acrylic -glass structure. For the three-layer acrylic-glass- 
acrylic structure, the electric fields are mostly shielded by the 25 
glass plate. Therefore, adding an acrylic plate beneath the 
glass plate does not result in a significant change in the sensor 
output signal. 

4.2 Detection of a Localized Anomaly in a 30 

Three-Layer Structure 

Water intrusion has been a persistent problem for compos- 
ite structures on aircraft. The freezing and thawing of intruded 
water in radomes and honeycomb sandwich flight controls 35 
can lead to disbond and structural failures. For this reason, 
water ingression experiments based on a sandwich structure 
were conducted to demonstrate the sensor’s capability of 
detecting water intrusion in radome structures. The sandwich 
panel used in the following water ingression tests, shown in 40 
FIG. 11 , has a paper and resin honeycomb core covered with 
fiberglass skins and closely resembles a real radome struc- 
ture. Table 1 gives the detailed properties of the sandwich 
panel. 

FIG. 1 1 shows the configuration for the coplanar capacitive 45 
sensor inspecting for 1 cc of injected water (4 honeycomb 
cells). The sensor scans from right to left on the test-piece 
surface, and the sensor output signal is read from the LCR 
meter screen. The solid line in FIG. 12 shows the sensor 
output signal for the configuration shown in FIG. 1 1 . It is seen 50 
from the solid curve in FIG. 12 that there are two peaks and a 
valley between them in the output signal. This phenomenon 
arises from the fact that the most sensitive region of the sensor 
is at the gap between its two electrodes. As the sensor scans 
over the water, the left gap of the sensor meets the water- 55 
injected area first. This results in a peak in the sensor output 
signal. As the sensor continues to move to the left and reaches 
the place where it is centered over the water-injected area, 
there is a decrease in the sensor output signal, due to the fact 
that the water is off the sensor’s most sensitive region. How- 60 
ever, as the sensor continues moving, its right gap then meets 
the water- injected area. As a result, there is another peak in the 
sensor output signal. When the sensor moves away from the 
water-injected area, its output signal returns to the baseline 
signal for the unflawed panel. 65 

In contrast, the dashed line in FIG. 12 shows the sensor 
output signal for 5 cc of injected water. In this case there is 
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only one peak in the sensor output signal, and the magnitude 
of the peak is approximately double that measured for 1 cc of 
injected water. This is due to the fact that the water-injected 
area in this case is laiger than in the previous case (20 hon- 
eycomb cells) . As the sensor scans from the right to the left, its 
left gap reaches the water-injected area first. Correspond- 
ingly, there is an increase in the output signal. As the sensor 
keeps moving, its left gap still lies over the water injected 
area, while its right gap starts to come into the water injected 
area as well. This leads to the maximum sensor output signal 
shown in FIG. 12. However, as the sensor continues moving, 
its left gap leaves the water-injected area first and the sensor 
output signal starts to decrease. When both gaps move out of 
the water-injected region, the sensor output signal returns to 
the baseline signal for the unflawed panel. 

5. Capacitive NDE System Using Sensor 

The present invention also provides for capacitive NDE 
systems which use the sensor. One example of such a system 
40 is shown in FIG. 13 . In FIG. 13 , a concentric coplanar 
capacitive sensor 10 is shown having a charged central disc 
14 . There is an outer annular ring 12 which is coplanar with 
and outer to the charged central disc 14 . In the embodiment 
shown, the rings 12, 14 are formed of copper and are placed 
on an insulative substrate, one example being a thin insulative 
film such as KAPTON® film. The sensor 10 is electrically 
connected to a capacitance measuring circuit 30 . The capaci- 
tance measuring circuit is electrically comiected to a proces- 
sor 32 which may be operatively connected to a display 34 . 
The system 40 may be housed in a housing 36 which may be 
handheld housing. 

In operation, the system 40 uses the concentric coplanar 
capacitive sensor 10 to determine a transcapacitance between 
a first electrode formed by the charged central disc 14 and a 
second electrode formed by the outer annular ring 12. The 
transcapacitance may then be used by the processor 32 as 
input to one or more models for a material under test. The 
transcapacitance sensed may be interpreted by one or more 
models for various purposes such as to determine permittivity 
ofindividual layers in a multi-layered structure orto allow for 
water detection (including water detection in radome struc- 
tures). The particular model used may depend upon the struc- 
ture being tested and its properties as well as the particular 
NDE testing being perfonned. 

Although a single probe system is shown, the present 
invention also contemplates the use of differential probes for 
optimal defect detection in the capacitive NDE system. 

Therefore, a concentric coplanar capacitive sensor and 
related methods and systems have been disclosed. The 
present invention contemplates numerous options, variations, 
and alternatives. For example, the present invention contem- 
plates variations in the materials used for the sensor, the 
specific size and geometry of the sensor, the type of structure 
being tested and the corresponding models for the structure 
under test, as well as other variations, options and alterna- 
tives. 

6. Rotationally Invariant Hand-Held Capacitive 
Probe 

A rotationally-invariant hand-held capacitive probe with 
concentric coplanar electrodes has been designed and built, 
FIG. 14 , motivated by defect detection in aircraft radome 
sandwich structures . Two versions of the probe, with different 
target penetration depths, have been tested. The sensors have 
the same outer diameter (25.4 mm) but different gap width 
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between the inner and outer electrodes. The probes were 
designed with the aid of the theoretical model previously 
discussed in which the capacitance is related to the electrode 
dimensions and the thickness and permittivity of each layer in 
a multi-layered dielectric test-piece. Experimental measure- 
ments of C with the probes in surface contact with one- and 
multi-layered dielectric test-pieces have been carried out and 
measured capacitance agrees with theoretical predictions to 
within 10%. The important parameters governing the pen- 
etration depth of this concentric capacitive sensor have been 
studied theoretically and experimentally by measurements on 
stepped Delrin® and HDPE slabs. Lift-off studies, both 
numerical and experimental, were carried out to investigate 
how lift-off affects measured C and the accuracy of the test- 
piece material properties when determined inversely from 
measured C. It is demonstrated that these hand-held probes 
are capable of detecting small embedded inhomogeneities in 
laminar structures, e.g., 1 cc of a low permittivity (low con- 
trast) injected fluid in a glassfiber-honeycomb-glassfiber 
sandwich panel that gives rise to AC-0.02 pF is clearly 
detected. On the other hand, significant impact damage in 
glass fiber composites was not clearly detected. These capaci- 
tive probes are especially promising for discontinuity detec- 
tion in sandwich structures. 

6.1 Theoretical Background 

FIG. 14A depicts a concentric capacitive sensor in surface 
contact with a five-layer dielectric halfspace. 

Hie capacitive sensor consists of an inner disc, radius s, 
and an outer annular ring, width t. The gap between these two 
electrodes is denoted g. A numerical model previously 
described provides a quantitative relationship between the 
sensor output signal, which can be measured, and the permit- 
tivity and thickness of each layer in the dielectric. In the 
model, the concentric sensor is considered to be infinitesi- 
mally thin while the test-piece is assumed to be laterally 
infinite. These assumptions are reasonable for electrodes that 
are relatively thin compared with the thickness of individual 
layers in the test-piece, and if the sensor is placed sufficiently 
far from the edges of the testpiece so that edge effects are 
negligible. In the theoretical model, the total charge Q on each 
sensor electrode is obtained from the calculated surface 
charge density and the sensor output capacitance computed 
from C=Q/V, where V is the potential difference between the 
two electrodes. For details of the calculation. 

6.2 Probe Assembly 

Two sets of concentric electrodes with different target pen- 
etration depths were fabricated by selectively etching a 
18-pm-thick copper cladding (14 mL standard) off a 25.4- 
pm-thick Kapton® film by photolithography (American 
Standard Circuits, Inc.). Both sets of electrodes have fixed 
outer diameter 25.4 nun (1 inch), which was selected as a 
workable dimension for a hand-held probe, but have different 
gaps and other dimensions as listed in Table 1. The charac- 
teristic capacitance listed in Table 1 is the calculated free- 
space capacitance for each sensor. The gap between the two 
electrodes and the width of the outer electrode are relatively 
small values and strongly affect the sensor output capaci- 
tance. In order to measure these values very accurately, a 
Nikon EPIPHOT 200 microscope was used that is capable of 
achieving precision of +5 pm for good calibration and 50x 
magnification. The sensor inner electrode radius was mea- 
sured using the “traveling microscope” method with accuracy 
+0.01 nun, due to its relative large dimension. It was found 
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that the fabricated dimensions are the same as the nominal 
values under such measurement accuracy. 


TABLE 1 


Dimensions and calculated free- so ace capacitance for sensors A and B. 


Inner 


Gap between 

Charac- 


electrode 


the 

teristic 


radius 

Outer electrode 

electrodes 

capacitance 


s (mm) 

width t (mm) 

g (mm) 

(pF) 

Sensor A 

10.67 ±0.01 

1.52 ±0.01 

0.518 ±0.009 

1.40 

Sensor B 

9.66 ±0.01 

1.527 ±0.008 

1.51 ±0.01 

0.99 


FIG. 15A shows the assembled capacitive probe, FIG. 15B 
shows the concentric electrodes, and FIG. 16 shows the com- 
ponents used to assemble the probe and components used in 
probe calibration. They consist of the following: a Rogers 
R04003® dielectric sensor substrate with thickness 
-, 0 0.31±0.01 mm, on which concentric electrodes are sup- 
ported; pins soldered to the electrodes; a BNC-to-receptacle 
adaptor that connects the pins to the BNC connector of an 
Agilent probe 16095 A; and an Agilent LCR meter E4980A 
that displays the measured capacitance. The entire sensor 
25 structure is enclosed in a two-part acrylic tube. Assembled 
parts A and B are shown in FIG. 15A with part B shown in 
detail in FIG. 15B. The acrylic tube was divided into two to 
facilitate calibration of the probe, i.e., removal of effects of 
the probe structure on measured capacitance. The two parts, 
30 which can be easily attached or detached, were connected 
together using plastic countersunk screws. 

6.3 Calibration Procedures 

35 .An effective calibration procedure removes the effect on 
the measured capacitance of all influences apart from the 
desired transcapacitance of the sensor. By comparing the 
probe measurement setup in FIG. 16 and the model used in 
numerical calculations (FIG. 14A), it can be seen that para- 
40 sitic capacitances that affect measurement results include: A) 
that from the cable connecting the LCR meter 72 to the BNC 
connector on the Agilent probe 70, B) that from the BNC-to- 
receptacle adaptor 68, C) that from the two receptacles 66 in 
which the two soldered pins 64 on the sensor 10 are inserted, 
45 and D) that from the two pins 64 themselves. The goal is to 
calibrate the whole system and take into account all the para- 
sitic capacitances up to the plane 62 shown in FIG. 16. In the 
LCR meter 72 measurement setup, the cable length option 
was set as 1 m. This setting automatically accounts for the 
50 parasitic capacitance due to the cable. In order to take into 
account parasitic capacitances from the BNC to receptacle 
adaptor 68 to the pins 64, open and short calibration steps are 
needed. Because the two pins 64 are soldered with sodder 60 
to the electrodes, as shown in FIG. 16 two identical pins 62 
55 were inserted into the ends of receptacles 66 during calibra- 
tion. Open and short calibrations were then performed on 
plane 62 according to the procedures provided in the LCR 
meter manual. All parasitic capacitances up to plane 62 are 
accounted for after calibration. However, effective permittiv- 
60 ity for the sensor substrate, as discussed below. 

6.4 Experiments on Laminar Structures 

As previously described, benchmark experiments measur- 
65 ing the transcapacitance of two concentric electrodes in con- 
tact with various large test-pieces showed agreement between 
experiment and theory of better than 4%. Similar experiments 
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are performed here to assess the level of agreement between 
theory and experiment for the hand-held probes, which is 
expected to be poorer due to the hardware associated with the 
hand-held probe that is not modeled explicitly. In order to 
account for effects from part B of the acrylic tube, an effective 
permittivity for layer 1, FIG. 14A, was introduced. This effec- 
tive permittivity was determined by placing the assembled 
probe in free space and measuring its capacitance. This mea- 
surement is considered in the numerical modeling as the case 
of a concentric capacitive sensor in surface contact with a 
one-layer dielectric (the sensor substrate) in free space. By 
assuming the thickness of layer 1 is the same as that of the 
sensor substrate and then varying its permittivity, a calculated 
probe output capacitance that agrees with the measured value 
to three significant figures was obtained. This permittivity 
value was subsequently assigned to be the effective permit- 
tivity of the sensor substrate with geometry shown in FIG. 
14B. The effective sensor substrate permittivity for sensor A 
was determined to be 3.47 while that for sensor B was deter- 
mined to be 3.31, at 1 MHz and room temperature. The 
effective permittivity values for both sensor configurations 
are greater than the substrate permittivity itself, 3.01 ±0.05, 
due to the existence of the acrylic tube part B (which has a 
relative permittivity of around 2.8). Hie effective permittivity 
of the sensor substrate for sensor A is greater than that for 
sensor B, because of the fact that sensor A has higher output 
capacitance values and influences from part B results in larger 
absolute changes in the capacitance for sensor A. Conse- 
quently, its effective substrate permittivity, which is inversely 
determined based on the output capacitance, is larger. These 
fitted sensor substrate effective permittivity values were sub- 
sequently used as inputs in the numerical model for the cal- 
culation of probe capacitances. 
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dielectric constant of the glass sheet was measured as 
5.62±0.05 and that of acrylic was 2.85±0.05. 

For all the measurements reported in this paper, the test- 
piece was supported 1 0 cm above a woodtop working table to 
5 approximate the free space assumption in the calculation 
model . The two probes with parameters listed in Table 1 were 
tested on five different laminar structures. The hand-held 
probes were pressed tightly against the test-piece surface to 
eliminate any air gap between the sensor substrate and the 
to test-piece. As can be seen from Table 2, experimental results 
agree with calculated results to within an average of 7% for 
sensor A and 9% for sensor B. Notice that absolute differ- 
ences in measured and calculated capacitance values for sen- 
sors A and B are similar in magnitude, and the greater relative 
15 differences observed for sensor B are due to the fact that its 
capacitance values are smaller. 

The agreement between theory and experiment of within 
10%, shown in Table 2, indicates that the structure of the 
probe give rise to some loss of quantitative accuracy, com- 
20 pared to the 4% agreement obtained in previous benchmark 
experiments for un-encased electrodes. Further, the calibra- 
tion process here is not perfect. For example, the electrical 
contact condition between the receptacles and the two sol- 
dered pins is not identical to that between the receptacles and 
25 the calibration pins. In addition, the soldered joints on the 
electrodes are not accounted for in the calibration process. 

6.5 Penetration Depth of Concentric Capacitive 
Sensors 
30 

In capacitive NDE, the penetration depth can be defined in 
terms of the sensor output capacitance [10, 14], Consider a 
concentric capacitive sensor in surface contact with a one 


TABLE 2 


Measured and calculated capacitance of hand-held probes 

in surface contact with various test pieces. 

Calculated C' (pF) Measured C' (pF) Relative Difference (%) 



Sensor A 

Sensor B 

Sensor A 

Sensor B 

Sensor A 

Sensor B 

One-layer acrylic slab 

2.75 

1.83 

2.58 ±0.01 

1.69 ±0.01 

-6.2 

-7.7 

One-layer glass slab 

3.57 

2.53 

3.26 ±0.01 

2.28 ±0.01 

-8.7 

-9.9 

Two layer glass over 
acrylic structure 

2.93 

2.03 

2.73 ±0.01 

1.87 ±0.01 

-6.8 

-7.9 

Two layer glass over 
acrylic-structure 

3.70 

2.63 

3.42 ±0.01 

2.37 ±0.01 

-7.6 

-10.2 

Three layer acrylic-glass- 
acrylic structure 

2.93 

2.03 

2.73 ±0.01 

1.87 ±0.01 

-6.8 

-7.9 

Average 





-7.2 

-8.7 


Relative differences are compared to the calculated capacitances. Uncertainty in measured C is 0.3% 


Measurements reported in this paper were performed at 
room temperature. The LCR meter operating frequency was 
set at 1 MHz so that the measurement error from the LCR 
meter was less than 0.3% for a 1 pF capacitance. At the same 
time, 1 MHz is low enough to be a good approximation for the 
electrostatic assumption made in the numerical model. 
Samples used in the benchmark experiments are one-, two- 
and three-layer test-pieces formed by combinations of acrylic 
and glass plates with lateral dimensions 30 cm by 30 cm. A 
digital thickness indicator with ± 1 pm accuracy was used to 
measure the plate thicknesses. The acrylic plates were 
2.39±0.02 nun thick and the glass plate was 3.02±0.01 mm 
thick. A Novocontrol Alpha Dielectric Spectrometer was 
used to provide an independent value of the dielectric con- 
stants of the samples at 1 MHz, as inputs to the model. The 


layer dielectric slab with permittivity £,. in free space (FIG. 
17A). The penetration depth D 10 of a concentric coplanar 
capacitive sensor is here defined by identifying the one-layer 
test-piece thickness T or which the capacitance is 10% 
smaller than its value when in contact with a similar but 
infinitely thick test-piece. When this condition is satisfied, the 
sensor penetration depth value D 10 is equal to the testpiece 
thickness T and is dependent on the permittivity of the test- 
piece. 

In other works, D 3 is defined as the penetration depth of 
capacitive sensors. Here we choose D 10 because the absolute 
difference in capacitance will be less than 0. 1 pF if the capaci- 
tance is less than 3 pF and D 3 is adopted, and such small 
changes in capacitance are hard to measure especially when 
noise is present. 


55 


60 


65 
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FIG. 17B shows the calculated sensor output capacitance 
as a function of the one-layer test-piece thickness and permit- 
tivity for sensors A and B and test -pieces with £,=2 and 5. It 
can be seen that, for a given test-piece permittivity, the sensor 
capacitance increases as the test-piece thickness increases 5 
and asymptotically approaches a constant value as the thick- 
ness becomes large. Further insight about the sensor penetra- 
tion depth is provided in FIG. 17C, in which the vertical axis 
is defined as 


k — Cool 

Difference % = X 100, 


C is sensor capacitance for a particular test -piece slab and 
is that as the slab thickness tends to infinity. Notice that the 
sensor output capacitance approaches C„ at different rates 
depending on sensor configuration and test-piece permittiv- 
ity. For a given test-piece permittivity, sensor B always has 20 
larger penetration depth than sensor A, because of its wider 
inter-electrode spacing. This agrees with our intuition. It is 
also shown that, for a given sensor configuration, the sensor 
penetration depth is larger for test-pieces with higher E, 
values. Test-pieces with larger permittivity values must have 25 
larger thickness Tin order to achieve a 1 0% difference in C for 
the slab with thickness T and a half space (see FIG. 17C), 
compared with test-pieces with smaller permittivity values. 

FIGS. 18A and 18B show measured capacitance as a func- , f) 
tion of test-piece thickness, for both stepped Delrin® 
(£,=3.82) and stepped HDPE (£,=2.65) slabs. FIGS. 18C 
and 18D show relative differences between the capacitance 
measured values and that on the 18-mm-thick test-piece, 
which approximates a half-space. It can be seen that D 10 of 35 
both sensors is greater for the Delrin® slab than that for the 
HDPE slab in accordance with the predictions of FIG. 19B, 
and both sensors’ sensitivity to test-piece thickness starts to 
decline as T increases. Fora given test-piece, D 10 for sensor B 
is greater than for sensor A again in accordance with predic- 40 
tions of FIG. 17C). Additionally, good agreement (to within 
an average of 10 %) between measured capacitances and 
numerical predictions is observed for test-pieces with permit- 
tivities and thicknesses in the range 2.65 to 3.82 and 1.50 mm 
to 3.12 mm, respectively. 45 

In summary, for a given sensor configuration, sensor pen- 
etration depth increases as test-piece permittivity increases. 

For a given test-piece material, sensors with wider inter- 
electrode spacing have higher penetration depths but smaller 
output capacitances. Therefore, a trade-off exists between 50 
sensor output signal and penetration depth. 

6.6 Capacitance as a Function of Probe Lift-Off 

How do lift-off variations affect the measured probe 55 
capacitance and the accuracy of test-piece permittivity values 
that may be derived from those measurements? The experi- 
mental arrangement for measuring C as a function of lift-off 
from the test-piece is shown in FIG. 15A. The test-piece was 
adjusted to be horizontal using a level. The lift-off between 60 
the hand-held probe and the test-piece was precisely con- 
trolled by pressing the probe tightly against the test-piece 
with fixed-thickness plastic shims acting as spacers in 
between. These plastic shims were then removed carefully, 
without moving the test-piece or the hand-held probe. This 65 
procedure helps to ensure that the plane of the electrodes and 
the test -piece surface are in parallel, avoiding probe tilt. The 
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thickness of the plastic shims was measured using a digital 
indicator and the resulting value considered to be the probe 
lift-off value. 

The capacitance of the hand-held probes as a function of 
probe lift-off was measured, and compared with numerical 
predictions. In the numerical calculations, the probe substrate 
was again assigned the effective value derived from measure- 
ment of the free space probe capacitance, and layer 2 in FIG. 
14A was assumed to be air with thickness equal to the lift-off 
value. The average difference between measured and calcu- 
lated values was 7%. The difference AC=IC ft »_ 0 «-C ai ,l is 
plotted in FIGS. 19A and 19B for measurements on PMMA 
and glass slabs respectively, whose parameters are described 
earlier in the section di scussing experiments on laminar struc- 
tures. Clift-off corresponds to the capacitance. 

Test-piece permittivity values can be determined inversely 
from measured capacitance values using the model. The mea- 
sured capacitances agree with the calculated ones the best 
when lift-off is large; since these situations are closest to the 
calibration environment of the probes. FIGS. 19C and 19D 
show the inversely determined permittivity values for the 
one-layered PMMA and glass, respectively. It is seen that 
when lift-off values are relatively small, the hand-held probes 
can characterize the test-piece material property fairly well. 
However, large lift-off values can result in inaccuracy in the 
inversely determined material permittivity information, even 
if the relative differences between the measured and calcu- 
lated capacitances are small. This is due to the fact that the 
hand-held probes are most sensitive to the region near the 
sensor substrate. When the lift-off is large a slight difference 
in measured capacitance can result in a large difference in the 
inversely determined test-piece permittivity. 

6.7 Detection of Embedded Inhomogeneities in 
Sandwich Structures 

In some structures, such as radomes, it is important that the 
electrical properties do not vary in an uncontrolled way. 
Imperfect repairs or damage followed by ingress of water can 
give rise to inhomogeneities with electrical properties that 
contrast with their surroundings. Here we investigate the abil- 
ity of the capacitive probes to resolve inhomogeneities of 
various size and permittivity embedded in a Delrin® plate and 
in the core of a glassfiber-honeycomb-glassfiber sandwich 
stmchire. 

Two rows of holes of different diameters, 2.5, 5.0, 7.5, and 
10 mm, were drilled in a 3.17-mm-thick Delrin® plate with 
permittivity £,=4.14. One row of holes was left empty while 
the other was filled with Paraffin wax (£,=2.1) to form zones 
with different permittivity contrasts. Table 3 lists the mea- 
sured capacitances when the contrast zones are positioned 
directly beneath the sensor gaps. It is seen that both probes 
were able to detect the air-filled holes of all sizes, whereas 
both failed to detect the 2.5-mm-diameter wax filled holes, 
due to the lower permittivity contrast between Delrin® and 
wax (around 2) compared with that between Delrin® and air 
(around 4.1). On the other hand, because of its deeper pen- 
etration depth, sensor B was found more capable of detecting 
embedded zones than sensor A; see relative differences in 
Table 3. 

In order to investigate the effectiveness of the hand-held 
probes in detecting inhomogeneities in sandwich structures, 
different amounts (1,3, and 5 cc) of water and olive oil (£,“8 
and 3 at 1 MHz and room temperature, respectively) were 
injected into the honeycomb core of a glassfiber-honeycomb- 
glassfiber structure (FIG. 20). Parameters of the sandwich 
panel are listed in Table 4. 1 cc of injected liquid corresponds 
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to 4 honeycomb cells with total surface area of 88 mm 2 , 
compared to the surface area of sensors A and B which is of 
507 mm 2 . FIGS. 21A to 21D show the capacitance measured 
as the hand-held probes scan a line directly over the cells 
containing the contrast agent. The measured probe signal 5 
strength is related to both the inhomogeneity permittivity and 
size. In particular, for the cases in which the injected liquid 
areas are smaller than the inner disc of the concentric sensor, 
two peaks in the output signal are observed for each measure- 
ment, due to the sensor gaps on each side of the sensor to 
responding to the inhomogeneity separately. On the other 
hand, when the injected liquid area is greater than the inner 
electrode size, a single peak in the measurement signal is 
observed due to both sides of the sensor being excited simul- 
taneously. As can be seen from FIGS. 21A-21D, the approxi- 15 
mate size of the inhomogeneity can be inferred from the shape 
of the measured signal and permittivity information can be 
extracted from the signal magnitude. 

In summary, the outstanding ability of the probes to detect 
low contrast zones smaller than the sensors themselves has 20 
been demonstrated. For example, both sensors detected suc- 
cessfully 1 cc of olive oil (G r “3) filling 4 cells in the honey- 
comb core of a laminar structure, indicating their potential 
application in defect detection in aircraft radome sandwich 
structures. 
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aged areas was 2.53±0.03 and for the impact damaged area 
was 2.54±0.01. This one example suggests that capacitive 
NDE is not suitable for characterizing impact damage in glass 
fiber composite. 

6.8 Conclusion 

Two hand-held capacitive probes with different target pen- 
etration depths have been built and tested. Following a cali- 
bration procedure that accounts for stray capacitances and the 
presence of the probe casing, which is not accounted for 
explicitly in the accompanying model, agreement to within 
1 0% between measured and calculated capacitances has been 
demonstrated for experiments on laminar structures. The pen- 
etration depth of concentric capacitive sensors lias been 
defined and studied both numerically and experimentally. For 
a given electrode configuration, the sensor penetration depth 
increases as test-piece permittivity increases. For a given 
test-piece, sensors with wider electrode spacing have larger 
penetration depths but lower capacitance values. The hand- 
held probes’ sensitivity to lift-off variations has been assessed 
numerically and experimentally. In order to acquire inversely 
determined material permittivities close to the actual values, 
small lift-off values are desirable because such measurement 
setups give rise to the best signal strength. This suggests that, 


TABLE 3 


Measured capacitance of hand-held probes on a Delrin ® slab with contrast zones. 


Hole 

diameter 

(mm) 

Measured capacitance 
on air filled holes foF) 

Relative diff. for air 
filled holes (%') 

Measured capacitance 
on wax filled holes (dF) 

Relative diff. for wax 
filled holes (%) 

Sensor A 

Sensor B 

Sensor A 

Sensor B 

Sensor A 

Sensor B 

Sensor A 

Sensor B 

2.5 

2.87 

1.91 

0.3 

1.0 

2.88 

1.93 

0.0 

0.0 

5.0 

2.85 

1.89 

1.0 

2.1 

2.87 

1.92 

0.3 

0.5 

7.5 

2.83 

1.86 

1.7 

3.6 

2.85 

1.90 

1.0 

1.6 

10 

2.79 

1.82 

3.1 

5.7 

2.82 

1.87 

2.1 

3.1 


The measured capacitances of intact areas are C = 2.88 pF and C = 1 .93 pF for sensors A and B, respectively. The relative difference is 
compared to the intact area capacitances. Uncertainty in measured C is 0.3%. 


TABLE 4 

Properties of the glassfiber-honecomb-glassfiber sandwich panel 


Parameter Value 


Core thickness 

7.62 

mm 

Skin thickness 

0.254 

mm 

Cell volume 

0.25 

cc 

Surface area of cell 

22 

mm 2 

Panel length and width 

298.45 

mm 


A glass fiber composite with dimensions 13.7 cmxl0.2 
cmx3 .24 mm was impact damaged on both sides by a dropped 
weight to generate a well-damaged area (of about 1 cm2 on 
each surface). Broken glass fibers and delamination were 
observed at the surface on both sides and assumed to exist 
throughout the whole thickness of the sample. Both sides of 
the sample where the weight was dropped were machined flat 
so that the signal of the capacitive sensor was from the inter- 
nal damage of the composite rather than from the surface 6Q 
indentation. 

Capacitive sensors A and B were used to assess the impact- 
damaged area. The capacitance values were then compared 
with those of the undamaged regions. For sensor A, the aver- 
age measured capacitance for undamaged areas was 65 
3.68±0.05 and for the impact damaged area was 3.67±0.02. 

For sensor B, the average measured capacitance for undam- 


if the probe is to be used for quantitative permittivity mea- 
surement. then calibration on a known test sample may be 
preferable to calibration in air. Experimental results show that 
the concentric capacitive sensors are unable to effectively 
characterize impact damage in glass fiber composites. The 
outstanding capability of the hand-held sensors in detecting 
relatively small contrast zones in one-layered and multi-lay- 
ered structures has been demonstrated experimentally, e.g., 1 
cc olive oil injection in glassfiber sandwich panel was suc- 
cessfully detected. 

The hand-held probes discussed here were built using 
readily available materials and components. In the future, 
some refinements can be made to the probe assembly in order 
to improve the agreement between measurement results and 
numerical calculations. For example, the probe test fixture 
and the BNC to receptacle adaptor can be replaced by a 
combined lead and sensor, thereby reducing parasitic capaci- 
tance. Additionally, the lead and sensor can be enclosed in a 
more compact rigid case that has fewer effects on the sensor 
signal. Thus, it should be understood that the present inven- 
tion contemplates numerous variations. In addition, the 
present invention contemplates variations in the materials 
used for the sensor, the specific size and geometry of the 
sensor, the type of structure being tested and the correspond- 
ing models for the structures under test. The present invention 
is not to be limited to the specific details of the embodiments 
described herein. 
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What is claimed is: 

1. A concentric coplanar capacitive sensor system, com- 
prising: 

30 a concentric coplanar capacitive sensor comprising (a) a 
charged central disc forming a first electrode, (b) an 
outer annular ring coplanar with and outerto the charged 
central disc, the outer annular ring forming a second 
electrode, (c) a gap between the charged central disc and 
35 the outer annular ring, and 

(d) the first electrode and the second electrode attached to 
an insulative film; 

a capacitance measuring circuit electrically connected to 
the concentric coplanar capacitive sensor for measuring 
40 transcapacitance between the first electrode and the sec- 
ond electrode foruse in evaluating a dielectric test piece; 
a processor operatively connected to the first electrode, the 
processor configured to use the transcapacitance as an 
input to a quantitative model of the dielectric test piece 
45 to determine inversely properties of the dielectric test 
piece; 

wherein the quantitative model provides a quantitative 
relationship between a transcapacitance measured with 
the concentric coplanar capacitive sensor and permittiv- 
50 ity and thickness of each of a plurality of layers in the 
dielectric test piece. 

2 . The concentric coplanar capacitive sensor system of 
claim 1 wherein the first electrode and the second electrode 
being formed of copper. 

55 3. The concentric coplanar capacitive sensor system of 

claim 1 wherein the charged central disc and outer annular 
ring provide rotational symmetry. 

4 . The concentric coplanar capacitive sensor system of 
claim 1 wherein the first electrode and the second electrode 

60 being electrically connected to a capacitance measuring cir- 
cuit for measuring transcapacitance between the first elec- 
trode and the second electrode. 

5 . A rotationally invariant hand-held capacitive probe, 
comprising the concentric coplanar capacitive sensor system 

65 of claim 1 . 

6 . The system of claim 1 wherein the dielectric test piece is 
a multi-layered planar dielectric structure. 
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30 


7. The system of claim 1 wherein the dielectric test piece is 
a multi-layered cylindrical dielectric structure with radius 
larger than the outer radius of the sensor by at least a factor of 
3 . 

8. The system of claim 1 wherein the dielectric test piece is 
a radome structure. 

9. The system of claim 1 wherein the dielectric test piece is 
a radome structure and wherein the processor is configured to 
use the transcapacitance to detect water or excessive inhomo- 
geneities caused by repairs in the radome structure. 

10. The system of claim 1 wherein the dielectric test piece 
is a radome structure and wherein the processor is configured 
to use the transcapacitance to detect inhomogeneities in the 
radome structure. 

11. The system of claim 10 wherein the inhomogeneities 
being caused by at least one of repair to the radome structure 
or impact damage. 

12 . The system of claim 1 further comprising a housing, the 
capacitance measuring circuit and the processor disposed 
within the housing. 

13. The system of claim 12 further comprising a display 
electrically connected to the processor, the display opera- 
tively connected to the housing. 

14. The system of claim 13 wherein the housing is a hand- 
held housing. 

15. A method of non-destructive evaluation, the method 
comprising: 

providing a concentric coplanar capacitive sensor; 

providing a quantitative model for the concentric coplanar 
capacitive sensor wherein the model provides a quanti- 


tative relationship between a transcapacitance measured 
with the concentric coplanar capacitive sensor and per- 
mittivity and thickness of each of a plurality of layers in 
a dielectric test piece; 

5 attaching the concentric coplanar capacitor sensor to the 
dielectric test piece; 

applying an input signal across the concentric coplanar 
capacitive sensor to produce an output signal; 

determining transcapacitance between the first electrode 
and the second electrode based on the output signal; and 

using the transcapacitance in the quantitative model to 
determine inversely properties of the dielectric test 
piece. 

! . 16. The method of claim 1 5 wherein the dielectric test piece 

comprises a plurality of dielectric layers. 

1 7 . The method of claim 1 5 wherein the dielectric test piece 
comprises a multiple layer aircraft radome structure. 

18. The method of claim 15 wherein the properties include 

20 a dielectric constant for each layer of a plurality of layers of 

the dielectric test piece. 

19. The method of claim 15 wherein the properties include 
thickness for each layer of a plurality of layers of the dielec- 
tric test piece. 

25 20. The method of claim 15 wherein the properties of the 

dielectric test piece include water in the dielectric test piece. 

21. The method of claim 15 where the properties include 
inhomogeneities in the dielectric test piece. 



